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ABSTRACT 

We present the results of C^^O(J=1-0) mapping observations of a 20' x 18' 
area in the Lynds 1204 molecular cloud associated with the Sharpless 2-140 (S140) 
H II region. The C^^O cube (a-(5-fLSR) data shows that there are three clumps 
with sizes of ~ 1 pc in the region. Two of them have peculiar red shifted velocity 
components at their edges, which can be interpreted as the results of the inter- 
action between the cloud and the Cepheus Bubble. From the C^^O cube data, 
the clumpf ind identified 123 C^^O cores, which have mean radius, velocity width 
in FWHM, and LTE mass of 0.36±0.07 pc, 0.37±0.09 km s-\ and 41±29 Mq, 
respectively. All the cores in S140 are most likely to be gravitationally bound 
by considering the uncertainty in the C^^O abundance. We derived a C^^O core 
mass function (CMF), which shows a power-law-like behavior above a turnover 
at 30 M0. The best-fit power-law index of —2.1 ± 0.2 is qu ite consistent with 



those of 



(120091) 



the IMF and the C^^ O CMF in the OMC-1 region by 



Ikeda fc Kitamura 



Kramer et al.l (Il998l ) estimated the power-law index of —1.65 in S140 
from the C^^0(J=2-1) data, which is inconsistent with this study. However, the 
C^^0(J=2-1) data are spatially limited to the central part of the cloud and are 
likely to be biased toward high-mass cores, leading to the fiatter CMF. Conse- 
quently, this study and our previous one strongly support that the power-law 
form of the IMF has been already determined at the density of < 10'^"^ cm~^, 
traced by the Ci*^O(J=l-0) line. 



Subject headings: ISM: clouds — ISM: molecules — ISM: individuals 140, L1204) — 
stars: formation 
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Introduction 



It is believed tliat new stars are born in dense molecular cloud cores. Therefore, the 
physical properties of the dense cores are thought to be closely related to the masses of 
new stars to form within them. Particularly, the dense core mass function (DCMF) has 
been considered as a key to understanding the stellar initial mass function (IMF). Many 
authors have investigated the DCMFs by using dense (> 10''"^ cm~^) gas tracers such as 
(sub)millimeter dust continuum emission, infrared visual extinction, and molecular line 
emissions having hi gh critical densitie s , for example, the H ^ ^CO^ (J=l-0) and N?.H^ 



(J=l-0) lines (e.g.. 



2007 



Motte et al. 



Rathborne et al. 



20091: 



1998 



Reid fc Wilson 



Ikeda et al 



2007 



2006 



Walsh et al. 



Nutter &: Ward-Thompson 



2007 



Enoch et al.l 



2008h. In 



nearby (< 500 pc) star forming regions such as Orion, Ophiuchus, the Pipe Nebula, 
Perseus, and Serpens, they found that the DCMFs seem to have power-law-like behaviors 
in high-mass parts as dN/dM oc M~'^, whose 7 values are very similar to that of the IMF. 
Based on these observational facts, they propose a hypothesis that the power-law form of 
the IMF has been already determined at the formation stage of the dense cores. 



Ikeda fc Kitamural (|2009[ ) have recently discovered a similarity between the power-law 
forms of the IMF and the core mass function (CMF) even in the tenuous (10^"*^ cm~^) cloud 
structures. They carried out mapping observations by using the C^^O (J=l-0) line, whose 



critical density is 2x10'^ cm ^ ( lYoshida et al. 



20101 ), toward the OMC-1 region in the Orion 



A cloud. They found the 7 value of the C O CM F of 2.3 - 2.4. T he 7 value is not only 



similar to those of the H^^C O+ DCMF (2.2 ±0.1; 



continuum DCMF (2.2±0.2; 



Ikeda et a. 



Nutter fc Ward- Thompson 



20071 ) and the 850 /im dust 



20071 ) within the un certainties, bu ; 



iMuench et al. 



also q uite consistent with that of the IMF of the Orion Nebula Cluster (2.2; 
2OO2I ). which is associated with the OMC-1 region. The agreement between the C^^O CMF 
and IMF 7 values suggests that, at least in the OMC-1 region, the power-law form of the 
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IMF with 7 > 2 has been maintained from the formation time of the tenuous structure 
with density of ~ 10^"^ cm~^. In addition to the Orion cloud, a representative of massive 
star formation, the similarity betwe e n the CMF and the IMF was reported in low-mass star 



forming regions by 



Tachihara et al. 



( 12002| ). The C^^O cores found in Taurus, Chamaeleon, 



Lupus, and other low-mass star forming regions show the CMF with 7 = 2.6 above 56 Mq. 



Kramer et al 



( 119981 ). however, found a significantly smaller 7 value of 1.7 for the S140 
and M17SW regions by using the C^^O (J=2-l) line with a high critical density comparable 



to that of (J=l-0). ,Wong et al.l (120081 ) also found 7 = 1.7 in the C^^^O (J=l-0) 



CMF of RCW 106. In this study, we focus on S 



distance to S140 of 910 pc (ICrampton fc Fisher 



dista nces to M17SW and RCW 106 are 2.2 kpc (IChini et al. 



40 fo r the following reasons. First, the 



1974r) is the near est one of the m; the 



19801 ) and 3.6 kpc (ILockman 



19791 ). respectively, and the two clouds are located in the Sagittarius arm. Since we discuss 
the power-law nature of the CMF on the basis of a comparison with that in OMC-1, it 
seems quite natural to first select the S140 region, located in the Local (Orion) Spur as 
the next step in our C MF study. Second, since the previous C^^O (J=2-l) observations 



of S140 (IJohnen 



19921 1 covered only the brightest 4' x 4' region around IRAS 22176+6303, 



we cannot exclude the possibility that their core identification was biased to high-mass 
cores, leading to the flatter CMF. Third, the higher transition of J=2-l with the transition 
energy of 15.8 K might prefer to pick up higher-mass protostellar cores, compared to the 
lowest transition of J=l-0. 

The aim of this paper is to re-estimate the power-law index of the CMF in the S140 
region, also identified as Lynds 1204, and to examine whether the similarity between the 
CMF and the IMF holds in the region or not. Our observations were done to cover the 



cloud as widely as possible, in order to avoid the possib 
C^^O (J=l-0) line emission, which is the same tracer as 



e spatial bias. We 



Ikeda fc Kitamural (120091 ) did, and 



emplo yed the 
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is suitable for deriving the ten uous CMF (e.g.. 



used th e clumpf ind algori t hm (jWilliams et al. 



Tachihara et al.ll2002l ). Furthermore, we 



19941 ). the same core identification method 



as that llkeda &: Kitamural (l2009l ) employed, to directly compare with their results. Note 
that the spatial resolution is high enough to evaluate the CMF at the distance to the S140 
region of 910 pc, as shown in ^ 

In ^ we describe the details of our observations. In ^ we show the overall spatial 
and velocity structures of the tenuous gas traced by the C^^O ( J=l-0) emission in the S140 
cloud, and identify three clumps with sizes of ~ 1 pc. In §l]we describe the identification 
of the cores in S140, and discuss their physical properties by comparing with those in 

the OMC-1 region. We show the C^^O ( J=l-0) CMF in ^ In this section, we compare our 
CMF with the previous J=2-l CMF, and discuss the similarity between the C^^O (1-0) 
CMF and the IMF. In ^we summarize our results. 



2. Observations 



The C^^O ( J=l-0) mapping observations were carried out in the period from January 
to February 2010 by using the Nobeyama 45 m radio telescope. Our map covered a 20' x 18' 
area, whos e center was selected to be 22'^20°^18^ 63°22'12" (J2000). Using the On-The-Fly 



technique (ISawada et al. 



20081 ). we swept the area by raster scan with a scan speed of 



the telescope of 35 arcsec s ^. To reduce scanning effects, we scanned in both the RA 
and Dec directions. A t the frequency of the C^^O (J=l-0) emission (109.782182 GHz; 



Ungerechts et al. 



19971 ). the half power beam width, A^hpbw, and main beam efficiency. 



f], of the telescope were 14" and 0.4, respectively. By the standard chopper-wheel method, 
the receiver intensity was converted into the antenna temperature T^, corrected for the 
atmospheric attenuation. We used an off position of 22^18^34*^, 63°2'12", where we could 
not the C^®0 emission above 0.09 K in Tt. At the front end, we used the 25-BEam 
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Array Receiver System (BEARS) in double -sideband (DSB) m o de, which has 5x5 bea ms 



separated by 41". 1 in the plane of the sky (jSunada et al. 



2000 



Yamaguchi et al 



2000|) 



The 25 beams have beam-to-beam variations of about 10% in both beam efficiency and 
sideband ratio. To correct for the beam-to-beam gain variations, we calibrated the intensity 
scale of each beam by observing the W3 region with a 100 GHz SIS receiver (SlOO) with 
a single-sideband filter. At the back end, we used 25 sets of 1024 cha nnel autocorrelato rs 



(AGs), which have a velocity resolution of 0.104 km s ^ at 110 GHz (jSorai et al. 



2000|). 



Since the data dumping time of the AGs was 0.1 s, the spatial data sampling interval on 
the sky plane was 3.5". The spatial interval corresponds to 0.3A6'hpbW) and is small enough 
to satisfy the Nyquist theorem. The telescope pointing was checked every 1.5 hours by 
observing the SiO ("^=1, J=l-0; 43.122 GHz) maser source T-Gep. Since the pointing 
uncertainty of the telescope remains as small as a few arcseconds below a wind speed of 5 
m s^^, we rejected the data taken under the condition that the wind speed averaged over 
one minute exceeds 6 m s~^. Gonsequently, the pointing accuracy is better than 3" for the 
data to be analyzed. 

To construct an a-(5-fLSR data cube with spatial and velocity resolutions of 22" and 0.1 
km s~^ in full width half maximum (FWHM), we used a Gaussian function as a gridding 
convolution function (GGF) to integrate the spectra which were taken with the very high 
spatial sampling rate of 3.5". We adopted 17". as the size of the GGF in FWHM. The 
resultant effective spatial resolution of the cube, A6'efr becomes 22". 0, corresponding to 0.1 
pc at the distance to the S140 region. Since the DSB system noise temperature ranged from 
285 to 357 K with a mean value of 321 K during the observations, we have an rms noise 



level of the data cube of 0.12 K in T^. 
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3. C^^O ( J=l-0) map of the S140 region 

In this section we describe the overall spatial and velocity structures of the C^^O 
(J=l-0) emission in the S 140 region. In addition, we identify three distinct clumps with 
sizes of ~ 1 pc, which is likely to be the natal objects for clusters. 



3.1. Spatial and Velocity Structures of the C^^O ( J=l-0) emission 



Figure [T] shows the total integrated intensity map of the C^®0 (J=l-0) emission of the 
S140 region. The total mass of the cloud is estimated to be 6600 Mq (see §4. II for details). 
One can see three distinct peaks, referred to as head, tail and filamentary clumps in this 
paper. The main body of the cloud is known as a cometary one, and our map clearly 
shows the cometary shape; the C^^O emission is the brightest toward the head clump, 
which faces the Sh-2 140 H ll region and is associated with IRAS 22176+6303, indicated 
by the cross mark in Figure [H The head clump is also known as an active cluster-forming 
region where numero us young stellar objects have been identified in infrared wavelength 



(IMegeath et al. 



20041 ). On the other hand, the tail clump seems streaming toward the 
north-east. In contrast to the head clump, no active star formation is known within the tail 
clump, although no sensitive infrared observations have been published yet. Since our map 
covers a wide r region by a factor of three than those of the previous molecular line mapping 



observations (IRidge et al. 



2003 



Higuchi et al 



20091 ). we found another filamentary clump 



located on the south-east of the cometary cloud. Note that the elongation of the filamentary 
clump seems perpendicular to the head-tail direction of the cometary cloud. Although the 
filamentary clump c an be seen in a previo us large-scale survey with the spatial resolution of 



several arcminutes (lYonekura et al. 



19971 ). our map is the finest one resolving 0.1 pc-scale 



cores. 
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As shown in Figure |2l we calculate the intensity-weighted mean LSR velocity, i.e., the 
centroid velocity, and veloci ty width in FHWM o f the C^^O emission at each spatial grid 



point in Figure [H following 



Yoshidaetal 



(l2010f ). The centroid velocity map illustrates 
that the centroid velocities in the cometary cloud range from —8.5 to —6.0 km s^^, while 
the filamentary clump has a more blue-shifted velocity range from —10.0 to —8.0 km s~^. 
The mean velocity width all over the cometary cloud is 1.0 km s~^, and the largest velocity 
width of 2.5 km s~^ is found toward the IRAS source, which is probably related to the 
star formation activity in the head clump. On the other hand, the velocity width of the 
filamentary clump is as small as 0.6 km s~^. However, the velocity width map apparently 
shows a somewhat large velocity width of ~ 2.5 km s"^ at a,S = 22'^20™43'^, 63°20'0" in 
the filamentary clump. As indicated in the position-velocity diagram of Figure [3]A, this is 
because a component belonging to the cometary cloud with a mean LSR velocity of ~ —7.5 
km s~^ overlaps with the filamentary clump with a mean LSR velocity of ~ —9.0 km s^^ 
along the line of sight. Note that the velocity width of each components is as small as ~ 
0.5 km s^^ at this point. 

The centroid velocity map indicates that there are red-shifted components of i'lsr 
= _7 _g 1^)21 on the periphery of the cometary cloud. Among the red-shifted 
components, the most prominent one can be found in the southern edge of the head 
clump. The red-shifted component of —6 km s~^, for example, is more easily recognized 
in the position-velocity diagram of Figure |3}3. Since the S140 region faces the Cepheus 
Bubble and the cloud exhibits the cometary shape, the region is likely to interact with the 



expanding bubb. 



(Abraham et al. 



e driv en by stellar winds and supernova explosions of members of Cep 0B2 



20001 ). Therefore, the red-shifted components around the cometary cloud 
can be interpreted as the outer envelope of the cloud that has been pushed toward the far 
side along the line of sight by the expanding bubble. However, the velocity widths of the 
red-shifted components of ~ 0.5 km s~^ never exceed that of the main body of the cometary 
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cloud. Since the expanding velocity of the bubble of 10 km (lAbraham et al.ll2000l ) is 
an order of magnitude larger than the thermal motion of the cloud, it is expected that the 
interaction is associated with shocks, leading to large velocity width. The small velocity 
widths of the red-shifted components might imply that the turbulent motions excited by 
the interaction in the red-shifted components have been already dissipated with a crossing 
time as short as 7x10"^ yr (see also a discussion in §4.2p . 



3.2. Identification of Clumps for Cluster Formation 



Lada fc Ladal (120031 ) showed that the clumps with sizes of ~ 1 pc and masses of the 
orders of 10^ - 10^ Mq are the natal objects for the stellar clusters. The three clumps 
in S140, the head, tail and filamentary clumps can be considered as the cluster-forming 
clumps, though active cluster formation has been found only in the head clump. This is 
because their sizes are about 1 pc and the masses are of the order of 10^ Mq. Actually, 
the masses of the head, tail, and filamentary clumps are estimated to be 2200, 2600, and 
1800 M0, respectively, from the total integrated intensities with t^lsr = —9.9 4.8 km 



s ^ of the areas above the 2 a noise l eve 
shown in §4.ip . Previously, iRidge et al. 



in F igure P (the details of the mass estimation is 



fl2003h and 



Higuchi et al. 



(120091 ) mapped the head 



clu mp by using the C O (J=l-0) line, and the clump mass was estimated to be 1900 Mq 



by 



Higuchi et al.l (120091 ) . which is roughly consistent with ours. 



Our interpretation that the three clumps in S140 are the site of cluster formation is 
also supported by the following virial analysis. Assuming the clumps are spherical, we 
derived the virial masses of 1100, 1300, 330 M0 with the uncertainty of a factor of 3 for the 
head, tail, and filamentary clumps, respectively. Since the clump masses are twice or more 
larger than the virial masses, it is most likely that the clumps are gravitationally bound 
and have the potential for producing clusters. In addition, the whole cloud of S140 is also 



-lo- 



in gravitationally bound state. By using the mean velocity width of 1.95 km s~^ in FWHM 
all over the cloud and the projected extent of the cloud of 3.7 pc of the cloud, the virial 
mass of the cloud is estimated to be 1900 Mq, which is by a factor of three smaller than the 
cloud mass of 6600 M©. 



4. C^®0 core catalog 



4.1. Identification of the Cores and Derivation of Their Physical Properties 



Following the C O CMF stu dy in the OMC-1 reg ion (llkeda &: Kitamura 



applied the clumpf ind algorithm (jWilliams et al. 



20091 ) ■ we 



19941 ) to the C^^O three-dimensional 



(a-(5-fLSR) cube data. The algorithm can work well with reasonable parameters to identify 
cores or cluraps, th ough several authors pointed out some shortcomings of the clumpf ind. 



Pineda et al 



(|2009[ ) examined the behavi or of the algorithm by changing the threshold 



level from 3 to 20 a, a wider range than IWilliams et al. 



(Il994j ) examined, and found 
that the power-law index of the mass function sensitively depe nds on the threshold for 



higher thresholds of > 5 a. However, 



Ikeda fc Kitamural (120091 ) demonstrated the weak 



dependence of core properties and CMF on the thresho l d in t he reasonable range from 2 to 



5 a levels, which had been also shown in 



Pineda et al 



(120091). Therefore, we adopted the 



threshold level for the algorithm of 0.24 K, i.e., the 2 a noise level o f the c ube data, which 



is recommended for identifying the core structure by 



Williams et al. 



(|l9M) and falls in the 



robust and reasonable range derived by llkeda fc Kitamural fl2009i ). We adopted the grid 
spacing of the c ube d ata of 22". 0, equal to A^cfi, i-e., full-beam sampling. This is because 



Williams et al. 



(jl994j ) determined the optimal threshold of the 2 a level for th e full- 



Deam 



Ikeda et al. 



fl2007h to 



sampling case. We also used the additional criteria introduced in 
reject ambiguous or fake core candidates whose size and velocity width are smaller than the 
spatial and velocity resolutions, respectively. 
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We identified 123 cores and estimated the beam-deconvolved radius -Rcore, velocity 
width in FWHM corrected for the spectrometer resolution (if core? LTE mass M^ore, virial 
mass Mvir, and mean density n of the C^^O cores. Table [T] shows the physical properties 
of the C^^O cores i n S14 0. The definitions of these parameters are the same as those in 



Ikeda fc Kitamural (1200^. Here we briefiy summarize the parameters specific to this study. 



We adopted AOcs = 22". as described in §21 the antenna efficien cy y of 0.4, and the 



19821 ). 



fractional abundance of C^®0 relative to H2, Xqisq of 1.7 x 10 ^ (IFrerking et al. 
The opt ical depth of the C^^ O line in the S140 region has been estimated to be smaller 



than 1; 



Higuchi et al. 



(l2009l ) derived the upper limit of the optical depth of 0.5, from 
the intensity ratio of the C^*^0 (J=l-0) to C^''0 (J=l-0) lines. We assumed that the 
excitation temperature Tex is uniform over the S140 region and is equal to the rotation al 



temperature of 24 K in the NH3 (1, 1) and (2, 2) observations by 



Higuchi et al. 



2009h . 



In the head clump, the temperature of > 20 K is reasonable because the c 



Sh 2- 140 H II region and numerous young stellar objects have been found (iMegeath et al 



ump faces the 



20041 ) ■ as well as the OMC-1 cloud. On the other hand, the other two clumps have not 
been well studied compared to the head clump. In the filamentary clump, two IRAS 
point sources 22192+6302 and 22196+6302 are detected, but the nature of the sources is 
unknown. Toward the tail clump, no signature of star formation has been found. Although 
our assumption of the uniform temperature cannot be validated for the two clumps, we 
found that a low Tex value of 10 K does not seriously affect our discussion, and therefore we 
adopted the assumption of the uniform Tex in this study. 



4.2. Physical Properties of the C^^O Cores 



We discuss the physical properties of the C^^O cores on the basi s of comparison wit 
those of the cores in the OMC-1 region of the Orion A cloud (jlkeda fc Kitamura 



th 
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20091 ). To fairly compare the S140 C^^O cores with the OMC-1 C^^O ones, we smoothed the 
OMC-1 data by a Gaussian function with a FWHM size of 37". 3 so that the OMC-1 cloud 
is put at the distance of 910 pc to S140. Furthermore, we randomly picked up 25 % of the 
original spectra to make the cube data for OMC-1 so as to have the same signal to noise 
ratio as that of the S140 data. Finally we re-identified 44 C^^O cores from the smoothed 
cube data in OMC-1. 

The left panel of Figure H] shows that the -Rcore distribution of the C^^O cores in S140 
has a single peak at 0.37 pc, which is much larger than the peak radius of 0.26 pc in OMC-1: 
the mean value of 0.36±0.07 pc for i?<.ore in S140 is larger than that in OMC-1 of 0.27±0.05 
pc. The Kolmogorov-Smirnov (K-S) test applied to the -Rcore histograms demonstrates that 
the distributions are considerably different from each other with a significance level of 1 %. 
The mass M^ore of the C^®0 cores in the S140 region tends to be also larger than that in 
OMC-1, as shown in the middle panel of Figure IH the mean value in S140 of 41±29 Mq is 
twice larger than that in OMC-1 of 18±9 Mq. In addition, the K-S test shows that the two 
distributions are considerably different from each other. Although only one broad peak at 
~ 20 Mq seems to exist for OMC-1, we can see two broad peaks at ~ 20 Mq and ~ 50 M© 
for S140. 

The above comparisons suggest that the large (-Rcore > 0.4 pc) and massive (Mcore > 
50 Mq) cores exist only in the S140 region. Furthermore, the mass fraction of the cores 
in S140 is very large: the total mass of the C^'^0 cores is 5000 Mq, 75% of the total mass 
traced by the C^^ O emission (see ^3.21). T his fraction is twice larger than that of the C^^O 



cores in OMC-1 (Ikeda fc Kitamura 



60 % for the Hi=^CO+ cores f Ikeda et al. 



2009), but is ro ughly consistent with the fraction of 



20071 . 



20091 ). The large massive cores and the 



high mass fraction of the cores in the S140 region are likely to be caused by high column 
densities in the region. Actually, the C^^O ( J=l-0) peak intensity of 2.4 K in in S140 is 
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significantly larger than that in the OMC-1 region of 1.6 K. It is possible to interpret that 
the accumulation of interstellar matter occurred owing to the interaction between the S140 
cloud and the expanding Cepheus Bubble, as suggested by the cometary shape and velocity 
features of the cloud as described in §3.11 

The n range of 1.0 to 9.8x10'^ cm~^ for the C^^O cores in S140 is similar to that in 



OMC-1 of 2.3 to 9.3x10^ cm ^, and is consistent w ith the critical 



(J=l-0) line of ^2x10^ cm'^ flYoshida et al. 



2010) ■ In addition, 



density of the C ^^O 



Rathborne et al 



(120091 ) 



showed that the mean volume density of their C^^O cores in the Pipe nebula is 7.3x10'^ 
cm~^, which falls in our n ranges of (1.0-9.8) x 10^ cm~^. In spite of the large range of n, 
the traced the traced densities of < 10^ cm~^ is significantly smaller than those of the dense 
gas tracers used by the previous CMF studies mentioned in ^ 

In contrast to -Rcore and Mcore, the distribution of the velocity width dvcore in S 140 
is quite similar to that in OMC-1, as shown in the right panel of Figure |H The K-S test 
shows that there is no considerable difference with a significance level of 1 %. The mean 
value of dvcore in S 140 is 0.37±0.09 km s~^ and is almost the same as that in OMC-1 
of 0.38±0.12 km s~^. It is likely that the turbulent motions excited by the passing of 
the expanding bubble have been already damped, though the bubble could significantly 
increase the velocity width of the C^^O cores. Considering the total kinetic energy of the 
bubble is 2.7x10^° erg in H I ga s, the distance to the bubble center is 910 pc, and the 



radius of the bubble is 5 degree (lAbraham et al. 



20001 ). the kinetic energy injected into the 



C^^O cores with a mean radius of 0.36 pc can be estimated to be 1.7x10"^^ erg, leading to 
a large velocity width of 1.2 km s~^; such the large value of dvcore cannot be discerned in 
the right panel of Figure HI Here, we assume tha t the fraction of the injected kinetic ene rgy 



to be converted into turbulence is 0.01 to 0.05 (iMac Low fc Klessen 



20071 ). Since the expansion velocity of the bubble is 10.2 km s ^ (lAbraham et al 



2004 



Ikeda et al. 



20001). the 
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crossing time of the bubble over the cores is est imated to be as short as 7x10^ yr, much 



smaller than the age of 1.7 Myr for the bubble (j Abraham et al. 



20001) ■ Therefore, one can 



expect the dissipation of the turbulent motions within the cores at the present time. 



In the left panel of Figure [5] we show the dvcore-Rcore relation of the C^^O cores. Although 
the distribution of the S140 cores tends to shift toward larger radii, as shown in Figure SJ 
both the distributions in S140 and OMC-1 seem to be similar to each other. Actually, the 
best-fit power-law function for S140 is (dt;core/km s-^) = (0.65±0.07) x (i?corc/pc)°-^^^°-^^ 
with a small correlation coefficient of 0.39, and the best-fit one for the OMC-1 cores is 
(c?^^core/km s"^) = (0.55±0.21) X (i?core/pc)°'^^^°'^^ with a correlation coefficient of 0.22; the 
two best-fit functions are consistent with each other within the uncertainties. In the right 
panel of Figure [5] we show the correlation between the virial ratio, Myir/Mcore, and Mcore- 
Since the mean value of the virial ratio is 0.3±0.2 in S140, all the C^^O cores are likely 
to be gravitationally bound by considering the uncertainty in Xqisq of a factor 3. Both 
the data in S140 and OMC-1 show negative correlations; the best-fit power-law function 
of the S140 cores is (M^ir/Mcorc) = (1.27±0.25) x (Mcorc/M0)-°-^3±o.o6^ ^^^^i a correlation 
coefficient of 0.48. A similar negative correlation in the C^^O (J=l-0) data is also found 
in low-mass s tar forming regions of Taurus, Ophiuchus, Lupus, L1333, Chamaeleon and the 



Pipe Nebula ( iTachihara et al. 



20021 ). The power-law index of —0.43 i s shallower than th e 



value of —2/3 expected for "pressure-confined" structures derived by 



Bertoldi fc McKee 



(119921 ). This fact also suggests that the self-gravity is dominant in the S140 C O co res 



compared to the ambient pressure, unlike the p Ophiuchi case by 



Maruta et al. 



torn . 



5. C^^O CMF in the S140 region 



Figure [6] shows a CMF of the C^^O cores. The CMF has a turnover at around 30 Mq, 
and a power-law-like shape in the high-mass part above the turnover. Above 30 Mq, we 
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applied a single power-law function by considering the statistical uncertainties and found 
that the best-fit power-law index 7 is 2.1±0.2. However, the 7 value signifi cantly differs 



from that of the C^^O (J=2-l) CMF in the S140 region of 7 = 1.65±0.18 flKramer et al. 



19981 ). We consider that this difference is likely due to their limited observation area in core 



survey. As sho wn in Figure [T], it is apparent that the area of the C^^O (J=2-l) observations 



(IJohnen 



1992) covers only the central intense part of the head clump, where the C^^O 



intensities are relatively stronger than those in the other areas in the S140 region, and 
the core identification is likely to be biased to massive cores. To examine this possibility, 
we re-identifed 26 C^^O (J=l-0) cores within the C^^O (J=2-l) mapping area, and we 
re-estimated the CMF from the 26 cores, as shown in the left panel of Figure [71 Although 
our spatial resolution of 22'^ is tw ice coarser than that of the C^^O (J=2-l) study of 13" 



(IJohne 



1992 



Kramer et al. 



19981 ) and hence our statistical uncertainties are large in the 



low- mass part of < 30 Mq, the CMF shows a power-law-like behavior in the high-mass 
part. The bes t -fit 7 is 1.5±0.6 and is consistent with the previous value of 1.65 derived by 



Kramer et al. 



(119981 ). Furthermore, we estimated the CMF all over the head clump. There 
are 44 cores in the head clump and the best-fit 7 is 2.0±0.4, as shown in the right panel of 
Figure [71 Although the core numbers are not sufficiently large for the statistical analysis, 
the 7 value in the head clump tends to be larger than that in the central intense part of 
the head clump, and are consistent with that of all the C^^O cores in our observed region. 
Therefore, we conclude that the CMF should be estimated from the mapping data that 
covers at least one c lump, which is thought to be the natal structure of a stellar cluster 



(ILada fc Lada 



20031) 



Our 7 value of 2.1=b 0.2 is consistent with that of the C^^O CMF in the OMC-1 region 



(llkeda fc Kitamura 



20091 ) . This agreement confirms that our observations could resolve 
star-forming cores even at the large distance of ~ 1 kpc. A poor spatial resolution is one 
of the major causes of the underestimation of 7. To examine the dependence of 7 on the 



- 16 - 



spatial resolution, we created the smoothed OMC-1 C^^O data cubes by changing the 
effective resolutions, as described in §4.2^ and derived the 7 values for the smoothed cubes, 
as shown in Figure [3 It is clearly shown that our resolution of 22", corresponding to 0.097 
pc for S140, can correctly estimate the 7 value wi thin the uncerta ir ities. Actually, our 7 
value is consistent with that in the C^^O study by iTachihara et al.l ( 120021 ). having a spatial 
resolution enough to resolve 0.1 pc-scale cores. In addition. Figure [H] predicts that 7 is 
considerably underestimated for the case that A6cs becomes larger than 0.1 pc, which is 



the minimum radius of the cores in the OMC-1 regi on at the hig . 



Therefore, the small 7 value of ~ 1.7 in RCW 106 (IWong et al 



coarse spatial resolution of 0.78 pc at the distance of 3.6 kpc (ILockman 



lest resolution of 0.061 pc. 



20081) is hkel v due to the 



1979h . 



This study concludes that the power-law shape with 7 > 2 in CMF holds even in 



tenuous structures with the densi 
recent work by 



Ikeda fc Kitamura 



:i es of 10^ ^ cm ^ of the S140 region, in addition to our 
mQ% in OMC-1. Furthermore, the 7 values of the C^^O 



CMFs in S140 a nd OMC-1 are quite consistent with that of the Galactic field- averaged 



IMF of 2.3±0.7 f.Kroupa 



200 ll ). These observational facts lead us to the hypothesis that 
the power-law nature in the IMF originates in molecular cloud structures with densities of 
less than ~ 10^"^ cm~^. 

Our conclusion of the resemblance between the CMF and the IMF is consistent with 
the recent theoretical works showing that such a resemblance should be understood as a 
statistical relation, rathe r than a one - to-on e correspondence between a core and a star 
to be formed within it. 



Smith et al 



(I2OO9I ) examined the formation and evolution of 



gravitationally-bounded cores by SPH simulations and found that the initial masses of 
the cores have a poor correlation with the resultant stellar masses. However, the shape 
of the CMF is maintained throughout the simulations and the resultant IMF shape is 



quite consistent with the initial shape of the CMF (see also 



Chabrier &: Hennebelle 



2010h . 
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Furthermore, ISwift fc Williamsl (120081 ) presented a supporting result that the resultant 
IMFs are insensitive to the differences in theoretical models for star formation within cores 
and that the input shape of the CMF is kept until the IMF within the current observational 
uncertainties. 

Since our studies are limited to the two clouds of S140 and OMC-1 up to now, it is 
urgent to perform systematic studies of the tenuous CMFs in various star forming regions. 
Since the S140 and OMC-1 regions are massive star/cluster forming ones, it is interesting to 
explore the CMF shape in regions where no star formation activity is found because stellar 
feedback processes, such as outflow, stellar wind, and expanding H l l region may cause som e 



influence on the physical properties of the tenuous structures (e.g., iNakamura fc Lil 120071 ). 
Although the global S140 cloud and the individual C^^O cores are likely to be influenced 
by the Cepheus Bubble as shown in §3.11 and 14.21 we cannot examine wh ether the CMF 



Rathborne et al. 



shape is affected by the bubble or not only from our S140 data. Recently. 
( 120091 ) estimated the DCMF in the Pipe Nebula, where no active star formation occurs, and 
found that the DCMF shape including the slope in the high-mass part is quite consistent 
with that of the IMF derived in the Orion Nebula Cluster associated with the OMC-1 



cloud ( iMuench et al. 



2OO2I ). where the most active star formations occur in the solar 
neighborhood. This study suggests that the stellar feedback processes seem not to play a 
dominant role in determining the CMF shape. 

For future observations of the tenuous CMFs, we require the following two points on 
the basis of our results. First, mapping observations should not be spatially biased. We 
have shown that the lower value of 7 in the C^^O ( J=2-l) CMF previously estimated in the 
S140 region is likely to be caused by the limited mapping area only for the intense portion 
of the cloud. Second, observations should be done with a high spatial resolution enough to 
resolve the 0.1 pc-scale core structures, as shown in Figure [HI To further investigate the 
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CMFs over the Galactic disk, one should extend the studies in the solar neighborhood to 
those in regions where enviro nmental conditions, such as metallicity , inter stellar radiation. 



and turbulence, are different (IKrumholtz et al. 



2010 



Schmidt et al 



20101 ) ■ The first step 



would be toward the neighbor, Sagittarius/Perseus arms. Actually, in M17SW and ROW 
106 which are located in t he Sagittarius arm, smaller power-law indices of the CMFs of 



1.7 have been reported (IKramer et al. 



1998 



Wong et all 12008) . To do observations of 



cores in the Sagittarius arms, approximately 2.0 kpc apart from the sun, we need an angular 
resolution of 10" to achieve a linear spatial resolution of 0.1 pc. Furthermore, an extreme 
environment is known in the Galactic center region. The Arches cluster, associated with the 



Galactic center, shows the present-day mass function wi ; 



flatter than that in the solar neighborhood (jStolte et al. 



1 7 o f < 2, which is significantly 



20051 ). To resolve the 0.1 pc at 



the distance to the Galactic center of 8 kpc, we need a very high spatial resolution of 2". 
ALMA is one of the best facilities to allow us to carry out spatially unbiased observations 
with the spatial resolution as high as 0.1 pc in the distant Galactic regions and is expected 
to reveal the physical relationship between the CMF and the IMF in the Galaxy. 



6. Summary 

We have carried out the mapping observations toward the S140 region in the C^^O 
( J=l-0) line. From the C^^O map, which covers a 20' x 18' area in the region, we construct 
the C^^O core catalog. Our results and conclusions are summarized as follows. 

1. We identified three clumps in the C^^O map. The head, tail, and filamentary clumps 
are likely to correspond to the natal objects of the stellar clusters because the clumps 
are in gravitationally bound state and the masses of them is comparable to the typical 
one of cluster-forming clumps. We found red-shifted components at the periphery of 
the head and tail clumps, which can be interpreted as the pushed gas by the collision 
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of the Cepheus Bubble. 



2. We identified 123 cores from the C O cube data. The radius and mass of the 
C^^O cores in S140 tend to be larger than those of the C^^O cores in OMC-1 



(llkeda fc Kitamura 



20091 ). The differences suggest that the cometary cloud formed 



by the compression due to the collision of the expanding Cepheus Bubble. On the 
other hand, the velocity widths of the C^^O cores in the S140 and OMC-1 regions are 
fairly consistent with each other, possibly indicating that the turbulence excited by 
the bubble has already been dissipated. 

3. We have demonstrated that the power-law index in the high-mass part of CMF, 7, 
could be underestimated due to insufficient spatial resolution. Therefore, we should 
estimate the CMF with spatial resolutions smaller than the minimum sizes of cores. 
Furthermore, we have found that the 7 value could also be biased if a mapping region 
is limited to a part of a clump. 

4. The 7 values of the C^^O CMFs in the S140 and OMC-1 regions are q uite consistent 



20011). Therefore, it 



with that of the Galactic field-averaged IMF of 2.3±0.7 ( iKroupal 
is likely that the power-law nature in the IMF originates in molecular cloud structures 
with densities of less than ~ 10^"^ cm~^. 
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Fig. 1. — Total integrated intensity map of the C^^O (J=l-0) emission ( i^lsr= —9.9 — 4.8 
km s~^) in the S 140 region. The contour intervals are 0.25 K km s~^ (corresponding to 2cr) 
starting at 0.5 K km s^^. The gray scale bar is shown at the right-hand side of the panel. 



The white cross indicates the positio n of IRAS 22 



mapping region in C O (J=2-l) by 



76+6303. The white polygon shows the 



JohnenI (1l992l ). The dashed lines roughly indicate the 



boundaries of the head, tail, and filamentary clumps (see §3.ip . 
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Fig. 2. — [Left) Map of the intensity- weighted mean LSR velocity, i.e., the centroid velocity 
map for the C^^O (J=l-0) emission in the S 140 region. (Right) Map of the velocity width 
in FWHM. The color bar is shown at the top of each panel. In the left panel, the arrows 
indicate the cutting lines of the position-velocity diagrams in Figure [31 The cross marks 
on the arrows show the origins of the position axis for the position-velocity diagrams. The 
centroid velocity and the velocity width are calculated within the areas abov e the 3 a noise 



l evel: the thin solid polygons delineate the 3 a noise level in Figure [T] (see also 



Yoshida et al. 



feoiol)). 
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Fig. 3. — Position-velocity diagrams of the C^^O ( J=l-0) emission along the cutting lines A 
(left) and B (right) indicated in Figure [2l The origins of the position axes are indicated by 
the cross marks in Figure [2J The contour intervals are 0.25 K (corresponding to 2a) starting 
at 0.38 K (3 a). 
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Fig. 4. — Histograms of -Rcore (^e/t), Mcore [middle), and dvcore (right) of the C^^O cores. 
The open and gray histograms mean the C^'^O cores in the S 140 region and those in the 
smoothed OMC-1 data (see §4.2p . respectively. 



Table 1. Physical properties of the C^^O cores in the S140 cloud 
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Fig. 5. — Velocity width-radius (Left) and virial ratio-mass (Right) relations of the C^^O 
cores. The filled and open circles correspond to the S 140 and the smoothed OMC-1 data, 
respectively. The solid and dahsed lines show the best-fit power-law functions for the S140 
and OMC-1 data, respectively (see §4.2p . 
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Fig. 6. — C^^O {J— 1-0) core mass function in the S 140 region. The error bars show the 
statistical uncertainty of y/N, where N is the sample number in each mass bin. The solid 
line indicates the best-fit power- law function in the high- mass part above 30 Mq. 
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ig. 7 . — (Left) Mass function of the C^^O cores identified in the mapping region by 



Johnen 



19921 ) in the C^^O (J=2-l) hne. Note that the mass bins are the same as those in Figure 



[6l The sohd hne indicates the best-fit power-law function above 30 Mq, and the dashed hne 
a power -law function with th e index of —1.65, which was derived from the C^^O (J=2-l) 



data by iKramer et al. 



(Il998l ). {Right) Mass function of the C^^^O cores identified all over 
the head clump. The solid line shows the best-fit power-law function above 30 Mq, and the 
dashed one is the same as in the left panel. 
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Fig. 8. — Power-law index 7 of the C^^O (J=l-0) CMF, derived from the smoothed data 



cube in the OMC-1 region (jlkeda fc Kitamura 



20091), is shown function of the effective 



spatial resolution, A9es, in arcsec. The upper x-axis is in pc at the distance to the S140 
region. The error bars indicate the uncertainties in fitting a power-law function to the CMF 
i n the high-ma ss part. The horizontal dashed line shows the 7 value of the Salpeter IMF 



(ISalpetei 



I955I). 
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20 


50 


63 
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0.44 


0.30 


49.1 
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22 


20 


54 


63 


17 
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0.52 


0.45 
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0.44 


0.42 


72.0 


15.8 


0.2 


3.6 



Table 1 — Continued 



I.D. 




R.A. 






Decl. 






T* 

A, peak 




dVcore 


-^core 




Mvir/Mcore 


n 






.12000.0 






.12000.0 


km s i 


K 


aicscc 




km s i 


M . 


M . 




10-'' cm-'^ 


116 


22 


20 


57 


63 


27 


1.0 


-7.4 


0.50 


75.7 


0.33 


0.35 


12.3 


8.6 


0.7 


1.4 


117 


22 


21 
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63 


16 
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22 


21 


6 


63 


18 
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-9.2 
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0.46 


0.47 


99.1 
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0.2 
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22 


21 


16 


63 


14 
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18.0 


3.4 


0.2 


4.1 
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22 


21 


16 


63 


16 


1.0 
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1.72 
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0.33 
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62.8 
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0.1 


7.2 
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22 


21 


20 


63 


15 


39.0 
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80.2 


0.35 


0.36 


53.5 
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0.2 


5.1 
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21 


36 


63 


14 


11.0 
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1.19 


41.5 


0.18 


0.18 


5.1 


1.2 


0.2 


3.5 
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22 


21 


39 


63 


14 
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0.20 


9.9 


1.7 
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4.6 



Note. — The typical uncertainty of each quantity is as follows: For i?corei 0.04 pc, derived from the uncertainty in the estimation of the core 
projected area. For dvcore, 0.13 km s~^, corresponding to the velocity resolution. For Mcore, a factor of 3 (see text). For Mvir, a factor of 3, 
derived from the uncertainties in ilcore and dvcoie- For M^i^/Mcoie, a factor of 4, derived from the uncertainties in Mcore and Mvir- For n, a 
factor of 4, derived from the uncertainties in ilcore and Mcore. 



